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We probe the transition between superfluid and Bose glass phases using quantum 
quenches of disorder in an ultracold atomic lattice gas that realizes the disordered 
Bose-Hubbard model. Measurements of excitations generated by the quench ex¬ 
hibit threshold behavior in the disorder strength indicative of a phase transition. 
Ab-initio quantum Monte Carlo simulations confirm that the appearance of ex¬ 
citations coincides with the equilibrium superfluid-Bose-glass phase boundary at 
different lattice potential depths. By varying the quench time, we demonstrate 
the disappearance of an adiabatic timescale compared with microscopic parame¬ 
ters in the BG regime. 

The dynamics resulting from tuning, or quenching, quantum matter across a phase transition 
provide fundamental insights into the nature of many-particle systems (/). One approach to un¬ 
derstanding this problem is the Kibble-Zurek (KZ) scenario (2-4), which links the dynamical gen¬ 
eration of excitations as equilibrium is disrupted during a quench to the critical exponents of the 
equilibrium phase transition. Despite the prevalence of disorder in quantum matter, little is known 
about how disorder influences non-equilibrium dynamics and affects the KZ paradigm in closed 
quantum systems (5). Understanding dynamic and non-equilibrium properties of disordered quan- 
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turn materials is of paramount importance to applications such as quantum annealing of disordered 
spin systems to benchmark adiabatic quantum computing ( 6 ). In classical systems, disorder can 
have a profound impact and lead to the formation of glasses, which display quench phenomena 
such as aging, experimentally inaccessible equilibration times, and relaxation dynamics that are 
unrelated to equilibrium states (7). Whether or not analogous behavior occur in strongly interact¬ 
ing disordered quantum systems is an open question. 

Here, we use a quantum quench of disorder in an ultracold lattice gas to probe the superfluid- 
Bose-glass (SF-BG) quantum phase transition (Fig. 1). We show that the appearance of excita¬ 
tions generated by the quench is determined by the ground-state phase diagram via comparisons 
to quantum Monte Carlo (QMC) simulations. By varying the quench time, we demonstrate that 
this sensitivity to the many-particle ground state occurs despite the disappearance of an adiabatic 
timescale (compared with microscopic parameters) in the BG regime. To make these measure¬ 
ments, we create an atomic realization of the three-dimensional disordered Bose-Hubbard model 
(DBHM) using ultracold 87 Rb atoms trapped in a disordered optical lattice (8). The DBHM is 
a paradigm for strongly correlated and disordered bosonic systems, such as 4 He in disordered 
substrates like aerogels, disordered Josephson-junction arrays, and long-wavelength properties of 
superconducting electron pairs (9). In the DBHM, a strongly interacting SF undergoes a quan¬ 
tum phase transition into a BG when subjected to disorder. The BG phase exhibits the peculiar 
property of lacking long-range order while possessing infinite superfluid susceptibility (9), and it 
is therefore viewed as a gapless insulator with finite compressibility that arises from the presence 
of quasi-condensates, or SF puddles, embedded in an insulating background. Disordered ultracold 
atom gases have been used to indirectly measure the SF-BG transition via transport and coherence 
measurements in ID ( 10 ) and 3D ( 11 ) disordered lattices and in ID quasi-periodic lattices ( 12 , 13 ). 

In our experiment, we prepare a gas consisting of (27 ± 2) x 10 3 atoms cooled to a few 
nanoKelvin in a harmonic optical dipole trap. A disordered cubic optical lattice formed from 
pairs of counter-propagating A = 812 nm laser beams and a 532 nm optical speckle field is su¬ 
perimposed on the gas ( 14 ). The atoms experience a potential energy shift proportional to the 
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Figure 1: (a) Schematic phase diagram of the DBHM and quench. The gas is quenched from the 
BG to the SF regime by rapidly reducing (green arrow) the disorder strength A from A 0 to zero at 
fixed U/t, which is determined by the lattice potential depth s. Equilibrium configurations and the 
disordered lattice potential (false color) are shown at three values of A. For sufficiently high A, BG 
(blue) and SF (light gray) phases coexist in the trap, (b) Time sequence for the measurement. The 
lattice potential depth and disorder strength are shown using red and green lines, (c) Equilibration 
is disrupted during the quench and excitations are produced, which are measured in TOF images 
(grayscale). Images are shown for A 0 ~ 0.5 E R (i) and A 0 = 0 (ii) at s = 12 E R . The white 
ellipse marked on the residual from a fit to the image marks the fitted TF radius. For sufficiently 
high disorder, excitations such as vortices are apparent (red arrow) after the quench, while smooth 
profiles are obtained at low A 0 . 

speckle intensity, which varies randomly in space, leading to disorder in the Hubbard tunneling 
t, interaction U, and site occupation e energies. The DBHM we realize is characterized by the 
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Hamiltonian 


H = - ^ ^ - !) + ^ ^ muj2r lnii ( 1 ) 

<ij> i i i 

where i and j index the lattice sites and () indicates that tunneling occurs only between adjacent 
sites. In Eq. 1, hi is the number of particles on site i, bi (b\) removes (adds) a particle from site 
i, m is the atomic mass, u is the geometric mean of the trap frequencies, r,. is the distance to the 
center of the trap, and // is the chemical potential. We measure all energies in terms of the recoil 
energy Er = h 2 / 2 niX 2 « 170nK kn. The distribution of the Hubbard parameters, which are 
broadened around the values for the uniform system, are precisely known (8,15). The strength of 
the disorder is characterized by the average potential energy A associated with the speckle, which 
is approximately equal to the standard deviation of the distribution of site occupation energies. 
The lattice potential depth s (which controls U/t) and A are independently adjusted by tuning the 
power of the lattice laser and 532 nm light. The range of s we sample in this work corresponds to 
a strongly correlated, quantum depleted SF for A = 0. 

We probe the BG-SF transition by measuring the amount of excitation produced by quenching 
A at fixed s. The disorder strength is linearly ramped from an initial value of A 0 to zero in 30 
ms (Fig lb), which is slow enough to avoid creating excitations solely via the time variation of the 
spatially inhomogeneous disorder potential (14). Based on general arguments regarding the phase 
diagram in untrapped systems, the BH phase will appear in the low-density edge of the gas for 
sufficiently high A 0 (9). For stronger disorder, the BG-SF boundary moves inward, encompassing 
more of the atoms. Excitations produced by the quench are measured using time-of-flight (TOF) 
imaging. By imaging after a long (50 ms) period of free expansion, vortices and other excitations 
are transformed into modulations of the density profile and the measured optical depth ( OD). 
These excitations are visible in the characteristic images shown in Figs, lc and 2b. For low A, 
the density profile after the quench and TOF is smooth, while for high A, features consistent with 
vortices are present. 
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To quantitatively characterize the amount of excitation present after the quench, we measure 
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where O l3 is the measured OD at the pixel indexed by i and j within a mask set by a smooth 
fitting function / that is the combination of a Thomas-Fermi (TF) profile and a gaussian (which is 
approximately the equilibrium SF distribution). This method was previously used to observe the 
quantum KZ effect by measuring excitations generated via a quench between MI and SF states in a 
“clean” lattice (16). Data for s — 11 Er and A « 0-1 Er are shown in Fig. 2a. It is apparent that 
excitations are not generated by the quench until a threshold disorder strength is crossed, above 
which x 2 increases approximately linearly with A 0 . Similar threshold behavior is observed for all 
s we sample in this work. 

This threshold behavior is generally associated with a phase transition and the quantum KZ 
effect (2-4). In the KZ scenario, an adiabatic transition from a phase disordered (e.g., BG) to an 
ordered (e.g., SF) state is impossible because of diverging characteristic length and time scales. 
Dynamically traversing a quantum phase transition by tuning (or quenching) a Hamiltonian pa¬ 
rameter such as A necessarily leads to the formation of domains and excitations such as vortices 
that persist even after the transition is crossed (1). In this case, the excitations occur when the BG 
state is present in the gas and the SF-BG transition is crossed as the disorder strength is reduced 
during the quench. 

To connect the observed threshold disorder with the SF-BG transition, we carry out exact QMC 
simulations of the equilibrium system using the same trap and lattice parameters, atom number, and 
speckle disorder as in the experiment (14). For trap-free geometries in the thermodynamic limit, the 
BG is characterized by a vanishing superfluid order parameter and non-zero compressibility (17). 
In contrast, the trapped system we consider exhibits domains corresponding to SF and BG phases 
that we distinguish using the spatial extent of the condensate. The condensate is identified as the 
macroscopic occupation of a single-particle eigenstate that we can obtain from the single particle 
density matrix p x = YlijiHh) (18-21). 

For clean systems (A = 0) and U/t < 29.34 ± 0.02 at ultra low temperatures (ksT/Ylt <C 1), 
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Figure 2: (a)-(b) Results of quench measurements, (a) The observed y 2 as A 0 is varied at s = 
11 En and the piecewise linear fit (red line) used to extract the threshold disorder A t h are shown. 
The error bars show the standard error in the mean for the 6-12 measurements averaged at each A 0 . 
TOF images obtained after the quench are shown for A 0 = 0 (i), A 0 = 0.7 (ii), and A 0 = 0.95 E R 
(iii). (c)-(d) Results from QMC simulations, (c) The upper bound N BG /N on the BG fraction is 
shown as a function of A for s = 11 E R . The error bars show the standard error in the mean for 
the QMC statistical noise. The insets are three-dimensional contour plots of the highest (blue) and 
second highest (green) occupation eigenfunction of the single-particle density matrix for A = 0.05 
(upper left) and A = 1 E R (lower right), (d) Density slices through the trap center are shown for 
A = 0.1 (i), A = 0.6 (ii), and A = 1 E R (iii). The blue (green) regions are the SF (BG) domains, 
and the color bar shows the average number of particles on each site. 


a single condensate extends throughout the system that coincides with local superfluid density 
order parameter (14). As A is increased, this behavior changes and the extent of the macroscopic 


6 


























condensate shrinks, leaving behind regions devoid of coherence with it. Since the SF-BG transition 
is of the continuous type, phase coexistence is forbidden (22), and we identify these regions as BG. 
To illustrate this behavior, we show the two highest occupation eigenfunctions of p\ for s = 11 E R 
and A = 0.05 and 1 E R in Fig. 2c. At low A, all single-particle states are spatially overlapped with 
the SF domain, and the second highest occupied state results from interaction-induced quantum 
depletion. For sufficiently high A, however, this extended state is replaced by a spatially localized 
mode that corresponds to a non-macroscopic and locally coherent superfluid puddle characteristic 
of the BG phase. 

To compare with the measurements, we compute the BG fraction N bg /N as the fraction of 
atoms in regions without a macroscopic condensate present. This estimate is an upper bound at 
non-zero temperature because of thermal excitation. As shown in Fig. 2d, the BG as defined by 
this criterion emerges at the edge of the gas and grows in extent and number as A is increased. 
Typical behavior for N bg /N at s = 11 E R as A is varied is shown in Fig. 2c. Similarly to the 
amount of excitation created by the quench in the experiment, N/N bg is only non-zero above a 
threshold disorder, above which it increases approximately linearly with A. 

We construct the SF-BG phase diagram shown in Fig. 3 by estimating the threshold disor¬ 
der A th for generating excitations in the experiment and for BG to appear in QMC simulations 
using a piecewise-linear fit to data such as those shown in Fig 2. The fitting function assumes 
constant behavior for disorder strengths less than A t h and linearly increasing behavior charac¬ 
terized by the free parameters A th and a slope for disorder strengths greater than A th . Several 
important features of the phase diagram are evident. The threshold disorder A tb is weakly de¬ 
pendent on s, and the QMC and experimental results agree within the 40% systematic uncertainty 
in A 0 ; there are additional systematic and statistical uncertainties arising from finite temperature 
and disorder averaging (14). This agreement—which demonstrates that the quench dynamics and 
production of excitations in this strongly disordered system are sensitive to the ground-state, equi¬ 
librium phases—supports the quantum KZ scenario. Furthermore, the observed threshold behavior 
cannot be explained by mean field theory, which predicts that a BG appears for infinitesimal disor- 
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der (23). Finally, the decrease in A th at higher s (i.e., larger U/t), which cannot be accounted for 
by general classical percolation mechanisms, implies that interactions facilitate the transition from 
SF to BG. 



Figure 3: Phase boundary between SF and BG regimes. Every point is the result from a piece-wise 
linear fit to data at fixed s (such as those shown in Figs. 2a and c), and the error bars show the 
fit uncertainty. The black squares are the experimentally determined values of A t h from quench 
measurements, and the red circles are the QMC simulation results. 

We explore the dynamical timescale of the SF-BG transition by varying the quench time. In 
the KZ scenario for clean systems, the amount of excitation and heat produced during a quantum 
quench typically display power-law dependence on the quench time (4). The knowledge of how 
this changes in disordered systems is limited to one-dimensional spin chains, which show logarith¬ 
mic behavior (24, 25). In our experiment, the dependence of x 2 on the quench time r q when the 
SF-BG boundary is crossed is too weak to detect. Typical data are shown in Fig. 4 at s = 10 E R 
for the SF domain at a disorder strength just below A t h (Ao = 0.35 Er) and for the BG regime 
(A 0 = 1 Er). To avoid complications from decay of excitations during the quench, we determine 
the amount of excitation by measuring the temperature T of the gas after allowing rethermaliza- 
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tion in the trap for 150 ms (14). We show the fractional deviation in the temperature (T — T 0 ) /T 0 
in Fig. 4, where T 0 is the temperature of the gas without disorder applied, in order to normalize 
heating from the lattice laser light. For the SF domain, a characteristic timescale for excitations 
to occur is evident. The data fit well to a decaying exponential function with a time constant 
4.4 ± 1.6 ms. Quenches with r q much longer than this time do not produce excitation. In contrast, 
in the BG regime, an adiabatic timescale is absent for r q up to 140 ms, which is approximately an 
order of magnitude longer than the tunneling time h/t « 15 ms: the slowest microscopic timescale 
present in the DBHM Hamiltonian. We cannot explore larger r q because heating from the lattice 
light results in a loss of the signal-to-noise ratio. 



Figure 4: Fractional change in temperature of the gas after quenches with different quench times 
r q . The filled circles are for the SF phase with A 0 = 0.35 Er, which is smaller than the threshold 
disorder A t h- The open circles are for A 0 = 1 Er, which is greater than A^ and corresponds to 
the BG regime. The error bars shown the standard error in the mean for the 8-15 measurements 
averaged for each point. The solid line is a fit to an exponential decay. 


The emergence of an equilibration time much longer than microscopic timescales is reminis¬ 
cent of glassy behavior in disordered classical systems (7). Our understanding of dynamics in 
the DBHM for two and three dimensions is limited (2(5), since direct simulation is intractable for 
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experimentally relevant numbers of particles. Whether the long timescale we observe in the BG 
phase is associated with critical phenomena or if is is connecting solely with disorder-induced 
glassiness is thus an open question. More work, such as measurements of how correlations change 
during the quench (27), is needed to clarify how the KZ scenario is altered by disorder. 
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